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Qh Abstract. The renormalization group framework can be applied to Quantum Field Theory 

Q on curved space-time, but there is no proof whether the beta-function of the gravitational 

^ coupling indeed goes to zero in the far infrared or not. In a recent paper [1] we have shown 

^ that the amount of dark matter inside spiral galaxies may be negligible if a small running 

'— ' of the General Relativity coupling G is present {SG/Gq < 10~^ across a galaxy). Here 

^r^ we extend the proposed model to elliptical galaxies and present a detailed analysis on the 

K*' modeling of NGC 4494 (an ordinary elliptical) and NGC 4374 (a giant elliptical). In order to 

2~: compare our results to a well known alternative model to the standard dark matter picture, 

pq we also evaluate NGC 4374 with MOND. In this galaxy MOND leads to a significative 

CN discrepancy with the observed velocity dispersion curve and has a significative tendency 

cn towards tangential anisotropy. On the other hand, the approach based on the renormalization 

^ group and general relativity (RGGR) could be applied with good results to these elliptical 

,__! galaxies and is compatible with lower mass-to-light ratios (of about the Kroupa IMF type) . 
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1 Introduction 

Currently there is a large body of data coming from cosmological and astrophysical obser- 
vations that is mostly consistent with the existence of dark matter. Such observations also 
suggest that the hypothesized particles that constitute dark matter have very small cross sec- 
tion and typically travel much slower than light. These lead to the cold dark matter (CDM) 
framework, which is one of the pillars of the current standard cosmological model ACDM. 

It is not only tempting, but mandatory to check if such dark matter particles exist 
(by detecting them in laboratory based experiments, for instance) and also to check if the 
gravitational effects that lead to the dark matter hypothesis could follow from a more detailed 
and complete approach to gravity. The effects of pure classical General Relativity at galaxies 
have been studied for a long time and, considering galaxy kinematics, the differences between 
General Relativity and Newtonian gravity are negligible (for a dispute on the latter, see 
however Ref. [2]). 

One of the subjects which attract significant attention currently is modifications of 
General Relativity in order to replace either dark matter or dark energy. When all the 
astrophysical and cosmological data are considered, there is currently no particular model in 
this class that was proved to be as successful as the ACDM model, however diverse nontrivial 
achievements were accomplished by such new approaches to gravity; whilst the ACDM model 
has its own difficulties and problems, particularly at the galactic scale[3-6]. Probably the 
most well known success of the modified gravity approach relies on galaxy kinematics, where 



MOND[7-9] is largely the most cited example and achieved success in some areas where the 
ACDM results are at least unclear [10, 11]. There are examples in the cosmological realm as 
well. In particular, considering the evolution of linear perturbations and its comparison to 
the CMB and the LSS data, results either close or identical to ACDM results can be found in 
different scenarios [12-15] . The Bullet Cluster [16] was some times cited as a definitive proof 
on the existence of dark matter, however this system can also be modeled from a modified 
gravity perspective [17, 18]. These phenomenological results motivate the use of gravitational 
theories that are different from General Relativity in its standard form. 

Independent on whether gravity should be or should not be quantized, we know that 
the matter fields should. It is well known that the renormalization group can be extended 
to quantum field theory (QFT) on curved space time (e.g., Refs. [19-21]). In particular, 
concerning the high energy (UV) behavior, there is hope that the running of G may converge 
to a non-Gaussian fixed point in accordance with the asymptotic safety approach [22, 23]. 
Our present concern is, however, not about the UV completeness, but with the behavior of 
G in the far infrared (IR) regime. 

In the far IR regime of quantum electrodynamics, one finds classical electromagnetism, 
and hence no renormalization group running of its coupling constant. This behavior is in 
accordance with the Appelquist-Carazzone decoupling [24] (see Ref.[25] for a recent deriva- 
tion). In the case of gravity (in the context of QFT in curved space time) the same effect 
of decoupling has been obtained for the higher derivative terms in the gravitational action 
[26, 27]. Currently, it remains unclear whether the Einstein-Hilbert action coupling param- 
eter behaves as a constant at the far IR or not. Since in pure theoretical grounds with 
no additional hypothesis it is hard to advance in this direction, this possibility of General 
Relativity deviation has been developed on different grounds a number of times before, e.g. 
[28-34]. 

In [1] we presented new results on the application of renormalization group corrections 
to General Relativity in the astrophysical domain. Previous attempts to apply this picture 
to galaxies have considered for simplicity point-like galaxies (e.g., [32, 33]). We extended 
previous considerations by identifying a proper renormalization group energy scale jjl and 
by evaluating the consequences considering the observational data of disk galaxies. We 
proposed the existence of a relation between jjl and the local value of the Newtonian potential 
(this relation was reinforced afterwards from a different approach, by using a scale setting 
formalism [35]). With this choice, the renormalization group-based approach (RGGR) was 
capable to mimic dark matter effects with great precision. Also, it is remarkable that this 
picture induces a very small variation on the gravitational coupling parameter G, namely a 
variation of about 10~^ of its value across a galaxy (depending on the matter distribution). 
We call our model RGGR, in reference to renormalization group effects in General Relativity. 

The main purpose of this work is to extend the analysis of [1] (see also [32, 36, 37], [38]) 
towards elliptical galaxies, testing the RGGR approach in this context. This extension is 
also of value to future RGGR applications, for instance in cluster of galaxies. Many of the 
elliptical galaxies behave as spherically symmetric stable systems which are mainly supported 
by velocity dispersions, in sharp contrast with the disk galaxies, which are axially symmetric 
and are mainly supported by rotation velocity. To this end, we deduce the additional effective 
mass introduced by RGGR, detail some aspects of its general behavior on elliptical galaxies 
and present a detailed fitting using recent observational data of the galaxies NGC 4374 
(giant elliptical) and NGC 4494 (ordinary elliptical). Moreover we use exactly the same 



data of NGC 4374 to evaluate MOND and compare its results with the RGGR ones.^ The 
numerical evaluation of the velocity dispersion (VD) curves and related procedures use a 
program made by us on Wolfram Mathematica. 

2 A brief review on RGGR 

The gravitational coupling parameter G may behave as a true constant in the far IR limit, 

leading to standard General Relativity in such limit. Nevertheless, in the context of QFT 

in cuved space time, there is no proof on that. According to Refs. [32, 36], a certain 

logarithmic running of G is a direct consequence of covariance and must hold in all loop orders. 

Hence the situation is as follows: either there is no new gravitational effect induced by the 

renormalization group in the far infrared, or there are such deviations and the gravitational 

coupling runs as 

Ar<-^ /\/f2 

/3^_i = /,^!f_ = 2z. ^^^^^^ = 2i/Go 1. (2.1) 



Equation (2.1) leads to the logarithmically varying G{^) function. 
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where /io is a reference scale introduced such that G(/io) = Gq. The constant Go is the 
gravitational constant as measured in the Solar System (actually, there is no need to be 
very precise on where G assumes the value of Go, due to the smallness of the variation of 
G). The dimensionless constant z/ is a phenomenological parameter which depends on the 
details of the quantum theory leading to eq. (2.2). Since we have no means to compute the 
latter from first principles, its value should be fixed from observations. Even a small v of 
about ~ 10~^ can lead to observational consequences at galactic scales. Note that the first 
possibility, namely of no new gravitational effects in the far infrared, corresponds to i^ = 0. 
The action for this model is simply the Einstein-Hilbert one in which G appears inside 
the integral, namely,^ 

c^ f R 



^RGGRb] = Y^ / ^-^V^d'^x. (2.3) 

In the above, G is an external scalar field, it satisfies (2.2). For a complete cosmological 
picture, A is necessary and it also runs covariantly with the RG flow of G [32, 34, 39]. In 
the above, the A term was not written since its role is expected to be negligible for the 
galaxy internal dynamics, and it will not be used through the main part of this paper. In 
the Appendix A we numerically solve the variation of A inside the galaxies that are analyzed 
in this paper (NGC 4494 and NGC 4374) and show that albeit the value of A can increase 
many times inside a galaxy, it is far from sufficient to lead to any significative observational 
effect. To our knowledge, this is the first time that the variation of A in galaxies was directly 
evaluated in the context of renormalization group effects. 

There is a simple procedure to map the solutions from the Einstein equations with 
the gravitational constant Gq into RGGR solutions. In this review, we will proceed to find 

^We do not do the MOND analysis for NGC 4494 since it was previously analyzed in Ref. [9] and this 
galaxy does not constitute a hard test for MOND. 
■^We use the ( — \- ++) space-time signature. 



RGGR solutions via a confornial transformation of the Einstein-Hilbert action, and to this 
end first we write 

G = Gq + 5G, (2.4) 

and we assume 5G/Gq <^ 1, which will be justified latter. Introducing the conformally related 

metric 

G 
{jfiu = -^9iMu, (2.5) 

the RGGR action can be written as 

SnGGB[9] = S^n[9] + 0{6G^), (2.6) 

where S^h is the Einstein-Hilbert action with Go as the gravitational constant. The above 
suggest that the RGGR solutions can be generated from the Einstein equations solutions via 
the conformal transformation (2.5). Indeed, within a good approximation, one can check 
that this relation persists when comparing the RGGR equations of motion to the Einstein 
equations even in the presence of matter [1]. 

In the context of galaxy kinematics, standard General Relativity gives essentially the 
same predictions of Newtonian gravity. In the weak field limit and for velocities much lower 
than that of light, the gravitational dynamics can be derived from the Newtonian potential, 
which is related to the metric by 

*o = - (l + ^) ■ (") 

Hence, using eq. (2.5), the effective RGGR potential $ is given by 

^ ^ c2(5G , , 

^ = ^N + iyyr- (2.8) 

An equivalent result can also be found from the geodesies of a test particle[l]. For weak 
gravitational fields ^n/c^ ^ 1 (with <I>jv = at spatial infinity), hence even if 5G/Gq <C 1 
eq. (2.8) can lead to a significant departure from Newtonian gravity. 

In order to derive a test particle acceleration, we have to specify the proper energy scale 
/i for the problem setting in question, which is a time-independent gravitational phenomena 
in the weak field limit. This is a recent area of exploration of the renormalization group 
application, where the usual procedures for high energy scattering of particles cannot be 
applied straightforwardly. Previously to [1] the selection of /i oc 1/r, where r is the distance 
from a massive point, was repeatedly used, e.g. [28-30, 32, 40]. This identification adds 
a constant velocity proportional to v to any rotation curve. Although it was pointed as an 
advantage due to the generation of "flat rotation curves" for galaxies, it introduced difficulties 
with the Tully-Fisher law [41], the Newtonian limit, and the behavior of the galaxy rotation 
curve close to the galactic center, since there the behavior is closer to the expected one 
without dark matter. In [1] we introduced a /i identification that seems better justified 
both from the theoretical and observational points of view. The characteristic weak-field 
gravitational energy scale does not comes from the geometric scaling 1/r, but should be 
found from the Newtonian potential $Ar, the latter is the field that characterizes gravity in 
such limit. Therefore, 

- = f{¥]- (2-9) 



If / would be a complicated function with dependence on diverse constants, that would lead 
to a theory with small (or null) prediction power. The simplest assumption, // oc ^at, leads to 
^ ex 1/r in the large r limit; which is unsatisfactory on observational grounds (bad Newtonian 
limit and correspondence to the Tully-Fisher law). One way to recover the Newtonian limit 
is to impose a suitable cut-off, but this rough procedure does not solves the Tully-Fisher 
issues [32]. Another one is to use [1] 






(2.10) 



where $o and a are constants. Apart from the condition ^q < 0, in order to guarantee 
5G/Gq <C 1, the precise value of $o is largely irrelevant for the dynamics, since <l>'(r) does 
not depends on $o- The relevant parameter is a, which will be commented below. The above 
energy scale setting (2.10) was recently re-obtained from a more fundamental perspective [35], 
where a renormalization group scale-setting formalism is employed. 

The parameter a is a phenomenological parameter that needs to depend on the mass of 
the system, and it must go to zero when the mass of the system goes to zero. This is necessary 
to have a good Newtonian limit. From the Tully-Fisher law, it is expected to increase 
monotonically with the increase of the mass of disk galaxies. In a recent paper, an upper 
bound on fa in the Solar System was derived [36]. In galaxy systems, z^alcaiaxy ~ 10""^, while 
for the Solar System, whose mass is about 10"^" of that of a galaxy, i^a|soiar system ^ 10^^^. It 
shows that a linear increase on a with the mass (ignoring possible dependences on the mass 
distribution) is sufficient to satisfy both the current upper bound from the Solar System and 
the results from galaxies. Actually, in Sec. 5.2 it is shown that a close-to-linear dependence 
on the mass can also be found for elliptical galaxies by using the fundamental plane. 

Once the ^ identification is set, it is straightforward to find the rotation velocity for a 
static gravitational system sustained by its centripetal acceleration [1], 

Contrary to Newtonian gravity, the value of the Newtonian potential at a given point does 
play a significant role in this approach. This sounds odd from the perspective of Newtonian 
gravity, but this is not so from the General Relativity viewpoint, since the latter has no 
free zero point of energy. In particular, the Schwarzschild solution is not invariant under a 
constant shift of the potential. 

Equation (2.11) was essential for the derivation of galaxy rotation curves. Since elliptical 
galaxies are mainly supported by velocity dispersions (VD), the main equation for galaxy 
kinematics in this case will not be eq.(2.11), but eq.(3.6) with the mass M{r) given by eq. 
(3.16). 

3 Mass modeling of elliptical galaxies 

While disk galaxies are extended gravitational systems mainly supported by rotation, ellipti- 
cal galaxies are mainly supported by velocity dispersions (VD). There is good evidence that 
the elliptical galaxies dealt in this paper (and many others) are close to spherical systems, 
and we will consider this approximation in this paper. 

For stationary spherical systems without rotation, the mean velocity at a small cell 
centered at position r is (v) (r) = 0, while the mean square velocity satisfies {y'^\ = (v^), 



where and <p refer to the angular components of spherical coordinates. For such systems, 
from the colissionless Boltzman equation one derives the following Jeans equation [42], 

i^(fo?) + V? = -f^. (3.1) 

a. or ^ ' r or 

In the above, (T'l{r) = (?;^) (r) — (vr) (r) = (^v'^.'j (r) is the VD radial component at the radius 
r, £ is the local luminosity density, /3(r) = 1 — (Tg(r)/cr^(r) is the local anisotropy and ^ is the 
potential of the total force per mass that acts in the cell centered at r (i.e., — V<I>(r) = i^). In 
the absence of dark matter, and considering Newtonian gravity, $ would be the gravitational 
potential generated by the stars alone (we consider only elliptical galaxies with negligible 
amount of gas). 

It is a well known theorem of Newtonian gravity that for spherical systems 

a^^a^ (3.2, 

or r^ 

where Go is the Newton's constant and M{r) the total mass inside the radius r. Therefore 
the effect of a spherical dark matter profile is to replace the total mass from the stellar one 
M^{r) to M{r) = M^{r) + Mdm{r). Hence, from the knowledge of M^r), Mdm{r), £{r) and 
/3(r), one can solve the Jeans eq. (3.1) to find the radial VD contribution induced by the 
stellar mass (o"^^,), the radial VD contribution from the dark matter (o"^^^), and the total 
radial VD, which satisfies 

al = al^ + aldm- (3-3) 

The above decomposition of CTr is a consequence of ^ being linear on the mass and the 
left hand side of eq. (3.1) being linear on a'^. The case of non- Newtonian gravity will be 
commented latter. 

Astrophysical observations cannot measure o"^ since the VD seen from the Solar System 
is not the radial, but the line of sight VD (which we will also call "projected VD"). From 
straightforward geometrical considerations, the observed projected VD can be computed from 
[42] 

The procedure described above to first evaluate o"^ to then find a^ is computationally 
demanding (specially for models with constant but free /3). A significant computational time 
improvement is achieved by bypassing the computation of o"^ [43] , 

o, , 2Go f°° /r\ £(r)M(r) , , , 

"'^"^-mj, ^■(ii)^V^*- <=*•") 

where 



K{u) ^ lu'P-^ 



l-^)^^^^--(i'^-^4)-K^'^4.^ 



(3.7) 



B{x, a, b) = Jq 1°" ^(1 — t* ^)dt is the incomplete beta function and T is the Gamma function. 

The above concludes a short review on how to relate mass distribution to observed 
velocity dispersions for Newtonian gravity. Nevertheless, almost all the above can be applied 
to non- Newtonian gravity as well, the single exception being the relation between $ and 
M{r). A useful procedure to generalize the picture above to other gravity theories, is to 
introduce the concept of total effective mass ME{r) as follows, 

^'(r)r'^ q(r)r'^ 
MEir) = -^ = ^^, (3.8) 

where g{r) is the norm of the local acceleration. 

As above discussed, the dark matter effect is to enhance the total mass inside a radius 
r from that computed from stars alone (M=|,(r)) to M{r) = M^,{r) + Mdm{^)- In the case 
of RGGR and MOND (without dark matter), a similar phenomenology is achieved since 
the potential $ is enhanced due to a change on the gravitational theory. The latter can be 
interpreted as the effect of an additional effective mass such that the total effective mass is 

ME{r) = M,ir) + M^acnir), (3.9) 

or, in the case of MOND, ME{r) = M^{r) + Mmond(?^)- 

For MOND, the acceleration g felt by a test particle due to the gravitational force 
of a spherical mass distribution is only equal to the Newtonian one for sufficiently large 
accelerations {g ^ ao). In this framework, the physical acceleration g is related to the 
Newtonian one {gn) by 

5/"(^)=5iV. (3.10) 

The function /i (MOND's interpolating function) is such that /u(x) = 1 for x ^ 1 and 
^(x) = X for X ^ 1 [7, 8]. This introduces a certain ambiguity to this framework, since each 
choice of the function /i leads to physically different results. Fortunately, for many sounding 
choices of //, and for the purposes of some physical tests, the differences are negligible. For 
concreteness, we will fix /x as one of the most cited interpolating functions, namely the one 
proposed in Ref. [44] which reads 

mW = ^. (3.11) 

With the above, the eq. (3.10) can be explicitly solved and it leads to the following additional 
effective mass contribution, 

MMOND(r) = ME{r) - M,{r) = ^M.(r) Jl + ^^ - ^. (3.12) 

The constant ag is assumed to be an universal constant for MOND (i.e., all galaxies should 
be subjected to the physical effects induced by the same value of ao). From a best fit to 
the kinematics of diverse disk galaxies, and considering the /i(x) function above, its value 
was found to be ao = 1.35 x 10~^cm/s^[45]. — If we had adopted the original /i function 
proposed in [7, 8], whose ao value is slightly lower according to disk galaxies data [45], a 
worse concordance with the NGC 4374 galaxy would be found. 



In the RGGR case, in order to find its effective additional mass, we start from the 
expression for the acceleration of a test particle in the weak field and low velocity limits, 
which was found in [1],'^ 

<GGa - ^'n (l - ^) , (3.13) 

where $Ar stands for the Newtonian potential of all the relevant matter. In the case of 
elliptical galaxies, assuming no dark matter, it is essentially the Newtonian potential from 
the stars (i.e., $7v = ^*)- 

Equation (3.13) is the analogous of the eqs. (3.10, 3.11), since it establishes the relation 
between the physical acceleration (considering its theoretical framework) to the Newtonian 
one. Contrary to MOND, eq. (3.13) poses no fundamental acceleration scale. In particular, 
there is no MOND's n{x) interpolating function that is capable to lead to the relation given 
in eq.(3.13). 

For spherical systems, 

Koo.ir) = ^^ (i + ^.^rr-^ — ] (3-14) 

2, 



GoM^r) c^av 



r 



' ^r+^^i-p.ir'ydr'- ^^'^^^ 



Hence, the effective additional mass of RGGR for spherical systems is given by 

2 2 

T OIVC T 

4 Observational data and numerical procedures 

4.1 Observational data 

In order to disclose the RGGR consequences to real elliptical galaxies, we selected two ellip- 
tical galaxies, an ordinary (NGC 4494) and a giant one (NGC 4374). The analysis of both 
of them are based on recent data, which in particular include extended line-of-sight velocity 
dispersion (VD) data obtained by both long-slit observations (inner radii) and by the VD of 
Planetary Nebulae (PNe). 

In Table 1 we list the main global properties of each galaxy which were relevant for the 
evaluation of their mass models. The expected stellar mass-to-light ratio is based on previous 
analysis that use the Kroupa IMF [50, 51], see also Ref. [52]. We have not added to the table 
above, but /3 estimates can be found in Ref. [52], where it is found that both of the galaxies 
are close to isotropic (—0.5 < j3 < 0.5). Converting to the V band the T^k's derived in Refs. 
[53, 54], it is fair to assume that, for NGC 4374, {^ *) ^,,^^^ ^y^^ ~ (4.5 ± l.O)M0/L0y (see 
also Ref. [48]). For NGC 4494, we use (T.,)^^^^^^ ^^^ ~'(3'^8 ± l.Q)MQ/LQy [49] (it should 
be noted that this galaxy kinematics is not compatible with the Salpeter IMF [52]). It is not 
our purpose to be ultimately precise on the Kroupa IMF, the essential issue is to verify that 
RGGR has a tendency of being in conformity with IMF's that lead to significantly lower T^,'s 



^In the presence of matter as a fluid this relation is not vaHd as an equahty, but within a good approximation 
(since the G variation is very smaU), the relation (3.13) holds. See Ref. [1] for further details. 
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17.1 
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113.5" 


6.11 
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[46-48] 


NGC 4494 


15.8 


2.64 


3.8 


48.2" 


3.30 


273" 


[47, 49] 



Table 1. 3) Luminosity in the V band. 4) The expected V-band stellar mass-to-light ratio considering a 
Kroupa IMF, see comments in text. 5) The effective radius. 6) The Sersic index. The last two values were 
found from Sersic profile fits to the observational surface brightness along the projected intermediate radius 
(for details see [48, 49]). 7) The radius at which the Sersic extension to the observational surface brightness 
is implemented here. 



than the Salpeter one. This tendency could thus provide another important physical test for 
our proposal once stelar population models become more precise. 

For each of the above galaxies, we use the observational surface brightness data up to 
a the largest radius where the observation is trustworthy, and then extend it with its Sersic 
profile (analogously to the case of disk galaxies, where the extension is implemented with 
a simple exponential profile). The observational surface brightness profile, which we call 
Mobs(-R)) is determined from observations up to a certain radius. It is currently known that 
the surface brightness of many galaxies can be well approximated by Sersic profiles [55] . From 
/iobs one can fit the parameters Re and n to find the corresponding Sersic surface brightness 
of each galaxy, which we call ^s(-R)- The latter can be used to generate an extended surface 
brightness profile //cxt, which is composed by ^obs from the center up to a given radius and 
^s from that point out to a larger radius. 

The intensity associated with the Sersic profile reads (see Ref. [56] for a review), 

/(i?) = /oe-''"te) , (4.1) 

where R is the projected radius, n is the Sersic index and Re is the effective radius (which is 
defined as the radius that encloses half of the total luminosity). Since the total luminosity 
is given by L = 2it J^ I{R)dR? , the constant /3„ must satisfy r(2n) = 27(2n, /3„), where F 
and 7 stand for the gamma function and the (generalized) incomplete gamma function (i.e., 
7(a,6) = Jq i"~^e~*dt). Solutions for j3n can be found numerically. 

Since most photometric observations are stated in unities of mag/arcsec^, the following 
relation between surface brightness ^(i?) and the intensity I{R) is useful, 



/i(i?) = -2.5 log 
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Lq VlSOx 60x60 
where Mq is the magnitude of the sun in the appropriate band. 



+ Mq, (4.2) 



4.2 The stellar contribution to the observed VD 

To find the stellar contribution to the total projected VD, we proceed as follows: 



i) Find the luminosity density, which comes from the inversion of eq. (3.5) [42], namely 

iir) = -I r ^ffi=dR. (4.3) 



TT Jr y/R- 



From the above, a vector {i{ri)} is built. This vector is numerically interpolated to generate 
a numerical function for i(r) (this step is necessary to improve computational performance). 
At this and similar procedures, the rj resolution is chosen such that the errors between the 
original function and the one build from the interpolation are no more than 1%. 

ii) Using a linear relation between stellar mass density and luminosity density, /9*(r) = 
T* i{r), it is straightforward to build the total stellar mass inside the radius r, M^{r). 

in) We use £, M* and the eq. (3.6) to build the numerical table {a'^{Ri, I3j)/T^}, with 
Pj G [—1, 1]. In general, f3 € (— oo, 1], but physical considerations (based on galaxy formation 
or reasonable restrictions on the distribution function) disfavor large negative (tangential) 
anisotropy, while favors either isotropy or radial anisotropy [48, 49]. 

The three steps above determines Up*, apart from T* and /3, which will be analyzed 
latter. 

4.3 The non-Newtonian contributions to the observed VD 

To find the RGGR additional contribution to the projected VD, which we call (Tprggr(-R) 
and satisfies 

we re-do the steps ii and Hi from above with Mrgqr (see eq. (3.16)) in place of M^. With 
this procedure, the total projected VD for RGGR is completely determined except for f3 and 
two constants that appear linearly inside a?,^ and o"prggr' namely: the stellar mass-to-light 
ratio T* and the constant I'a. 

The effective additional mass provided by MOND (eq. 3.12) has a simple form, but it 
depends on different powers of T*. Hence, in the step Hi described previously, the numerical 
function (Tp mono is found from the three dimensional interpolation of the following table, 
{(Tp^Qj^j3(T^,j, i?j,/3fc)}, where T,,, was constrained to be in the range'* 1 < T=k/(M0/Lq^) < 
10, while R and /3 have the same resolutions and range as in the RGGR model described 
above. 

Before proceeding, we remark that the stellar component modeling used for MOND is 
identical to the modeling used for RGGR and Newtonian gravity. In Ref.[57] it is argued in 
favor of the use of Jaffe profiles in the outer parts of the baryonic mass profiles of ellipticals. 
Here this route is not pursued. 

5 Elliptical galaxies within RGGR without dark matter: general aspects 

5.1 Sersic profiles 

For spherical stationary systems, the main differences between Newtonian gravity and the 
RGGR one are encoded in the expression for Mrggr (3.16). From a given matter density 

''No physical issues are expected due to this constrain for this galaxy in this band. 
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Figure 1. The RGGR additional mass contribution A/rggrj as a function of the deprojected radius 
r, for two Sersic profiles. One of the profiles has Sersic index n — 6.1 (solid line) and the other with 
n = 3.3 (dashed line). Many of the elliptical galaxies have intermediate Sersic indices. For av ~ 10" 
the additional mass Mr^ggr is about the same order of the baryonic mass in a galaxy. 
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Figure 2. The projected VD contribution from RGGR (ctprgghX^)) for a matter distribution given 
by the Sersic profile (4.1) with n — 6.1. The thickest (green) curve corresponds the isotropic case 
^ = 0. The others are associated to the following /3 values: -1.00, -0.75, - 0.50, - 0.25, 0.25, 0.50, 0.75, 
1.00 (from light gray to black, respectively). 



p{r), it is straightforward to derive the additional effective mass posed by RGGR. In the end, 
the single uncertainty relies on the constant a, on which further comments on its behavior 
will be presented in the next section. 

Since the matter density of diverse elliptical galaxies are fairly proportional to a depro- 
jected Sersic profile (which, for n = 4, coincides with the de Vaucouleurs profile), here we 
will consider the consequences of RGGR for matter profiles deduced from deprojected Sersic 
profiles (for all radii). Two cases will be explicitly evaluated: the case n = 3.3 and n = 6.1 
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Figure 3. The same of Fig. 2, but witli Sersic index n — 3.3. 



(the same Sersic indices of NGC 4374 and NGC 4494). Such values are useful to the purposes 
of this section, since many ellipticals have Sersic indices that lie between these values. 

The RGGR mass profile M^QQ^(r) is shown in Fig. 1. The contributions of such profiles 
to the projected VD (see eq. (4.4)), including different assumptions for the anisotropic 
parameter /3, are depicted in Figs. 2, 3. 

Apart from the a role, from the eq. (3.16) one sees that Mrggr is not sensible to either 
the constant Iq in eq. (4.1) or the value of any constant mass-to-light ratio; it only depends 
on the effective radius Re and the Sersic index n. 

5.2 The role of a and the fundamental plane 

Here we point and illustrate that future analysis on the relation between RGGR and the 
fundamental plane [58, 59] can provide significant constraints on the correlation between a 
and the mass distribution. It is beyond the scope of this work to do a solid analysis on 
the RGGR consequences for the fundamental plane interpretation (interpretation which is 
currently an open problem within any framework). In particular, the status of such task for 
MOND can be checked in Ref. [60] and references therein. 

The existence of a plane in the space ((/)g ,Re,crQ) (the fundamental plane) is usually 
interpreted as a consequence of the virial theorem (in the context of Newtonian gravity) , see 
Ref. [61] for a review. Also, the data scattering can be further reduced around the theoretical 
surface if the Sersic index n is included as a new axis, a "fundamental hyperplane" [62]. It 
is admissible to use the same relation in context of RGGR, except that one should replace 
the virial mass My by the corresponding effective mass Mev (3.9) (considering stationary 
and spherically symmetric systems). Note that the total effective mass is not finite, but such 
feature is not a novelty, since well known dark matter profiles do display the same feature. 
One may define, for instance, the virial radius at the radius at which the local effective density 
is about 200 times the cosmological critical density. The details of such definition will not 
be relevant to the approach below. 

The theoretical plane in the ((/)e , Re, cq) space, which is derived from Newtonian grav- 
ity, together with the virial theorem and with constant M/L, is tilted in regard to the 
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observed one. The tilt can be interpreted as an evidence of dark matter (within the Newto- 
nian gravity context), since the necessary M/L variation to match theory and observation 
is larger than the expected from stellar population variations (e.g., [63]). Hence, one way to 
explain the tilt of the fundamental plane is from the assumption that 

focL^(/)f, (5.1) 

where the values of A and B depends on the values of the parameters that determine the 
fundamental plane (from observations). For instance, using the values of Ref. [64] one finds 
A = 0.31 and B = 0.02 (see also Ref. [61]). Part of the above scaling relation would be due to 
systematic stellar mass-to-light ratios {M^f/L) variations. For simplicity, we will here assume 
no systematic variation of M^,/L, hence 

At the virial radius, it is fair to use M* <^ M^qq^. Combining the last comment with 
eqs. (5.1,5.2), 

^oc/o^+^i??V+^(n), (5.3) 

since (1)^ = Iq g{n)/(2TT) and L = lQR^g(n), where 






gin) = 27r^^^r(2n). (5.4) 



In Sec. 5.1 it is shown that M^QQ^{r)/a does not depend on /q, hence from eq. (5.3) a 
rough estimate on the variation of a with M^, can be found as 

aocMi+^+^/(i?e,n). (5.5) 

That is, for elliptical galaxies of about the same shape (i.e., with similar Re and n), a should 
increase with the stellar mass of the the galaxy faster than linear, but slower than quadratic. 
This is a quick procedure to evaluate the dependence of a with the galaxy parameters. The 
best procedure would be to use a large sample of elliptical galaxies to deduce that, and also 
unveil the function f{Re,n). 

6 NGC 4374 and NGC 4494 results 

6.1 Introduction 

Here a detailed numerical analysis of NGC 4374 and NGC 4494 is done. The main results 
are shown in Tables 2, 3, 4, 5 and Figs. 6,5, 8. The main purpose of this section is to show 
that recent data on elliptical galaxies is in conformity with the RGGR model, even assuming 
a negligible amount of dark matter. Moreover, we compare our results to MOND (which also 
constitute an original work presented in this paper). 

Many elliptical galaxies display a certain degree of rotation. Such rotation do not appear 
directly in the velocity dispersion (VD) data, but it also traces the galaxy effective mass. In 
order to mass model such galaxies without simply neglecting the rotational (sub-dominant) 
component, one can insert a compensation in the VD data. The data that is here called 
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Table 2. NGC 4374 results. (1) /3[o] indicates isotropic VD, /3[-i,i] indicates constant anisotropy with 
/3 £ [—1, 1], K.IMF is a reference to Kroupa IMF, and it means that the expected value of T* was used, see 
Table 1 and eq. (6.2). (4) and (5) are the chi squared and reduced chi squared parameters, for details see 
Sec. 6. 2. 
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Table 3. RGGR mass resuhs for NGC 4374. 



"line of sight VD" has aheady a compensation for rotation, whose computation was done 
in Refs.[48, 49]. We use the same data that these references use for mass modeling these 
galaxies. 

In the Newtonian gravity without dark matter framework, there are two parameters that 
the procedures described in Sec. 4. 2 do not fix, namely: the mass-to-light raio T* and the 
anisotropy parameter (3. There are theoretical expectations for T*, but they are considerably 
more uncertain than other astrophysical data. 

Contrary to MOND, the RGGR approach was not developed in order to remove the 
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Table 4. NGC 4494 results. See Table 2 for details. 
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K.IMF ± /3[o] 9.7 ±1.3 

K.IMF ± /3[_i,i] 10.0 ±1.6 

Table 5. RGGR mass results for NGC 4494. 



necessity of dark matter, but to explore a possible QFT effect on large scales that may have 
a significant impact on astrophysics and cosmology. In this paper, likewise in Ref. [1], we 
explore the extremum possibility of having no dark matter in galaxies and only standard 
gravitation with the running of G as given by eq. (2.2). It would be a leap of faith the 
assume that this single running of G could solve all the dark matter issues at all scales, and 
we do not expect that. On the other hand, the galactic picture we are describing can be more 
than a mere starting point, since it does work as reasonable approximation of a universe with 
a lower amount of dark matter, and of a warmer type. Considering this picture, dark matter 
itself would have a lesser role in galaxies, while the RGGR role would be the dominant one. 
Albeit the results here presented for NGC 4374 favor RGGR over MOND, a natural 
criticism would be that with RGGR one more parameter is being fitted than with MOND. 
The comparison done here is worthwhile since: i) had we systematically found for RGGR 
a worse observational concordance than for MOND, that would sign that RGGR had deep 
problems and that it should be dismissed; ii) in case MOND systematically deviates from 
the observational data (e.g., data from giant elliptical galaxies), it should be dismissed or 
modified, in) MOND and RGGR seem to display different tendencies for the stellar mass- 
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Figure 4. NGC 4374 mass models with Newtonian gravity (NG). In the larger plots at the left, the circles 
and stars with error bars refer to the observational VD, either from long-slit observations (the stars) or 
from planetary nebulae (the circles). The curves refer to mass models composed by a stellar component 
with Newtonian gravity (NG). The first model (i) assumes isotropy (/3 = 0), while the second (ii) assumes 
P G [—1, 1]. Models iii and iv use the expected mass-to-light ratio from the Kroupa IMF (K.IMF) as part of 
the data to be fitted (see eq. (6.2) and Table 1 for details). The four thin lines that appear in the lower left 
plot depict the theoretical error bars on the stellar velocity dispersion (the solid thin line refers to the model 
iii, while the dashed thin line to the model iv, which differ since they have different /3 values). The vertical 
dashed line signs the radius above which the observational data is considered for the fitting procedure (10 
arcsec) [48]. The four small plots at the right show the la and 2cr confidence levels for each of the models. 



to-light ratios, and thus developments on stellar population models may favor one over the 
other, in spite of any differences on the number of parameters. 

6.2 On the x^ minimization 

In order to find the best fit for each model, we proceed with a standard x^ minimization. 
We numerically search for the global minimum of 



2/"^ r^ nk \ S~^ I ^modc\{i^ *T^l:'<y2i ^i) CTobsi 



x'(T„ei,e2) = ^ 



i=l 



Ci 



(6.1) 



In the above, N is the total number of observational VD points, 0i and G2 stand for addi- 
tional parameters (e.g., /3 and av), a^bBi is the VD observational value at radius Ri, while Ci 
is its corresponding uncertainty. 
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Figure 5. NGC 4374 mass models with RGGR. The curves refer to mass models composed by the 
stellar component (inferred from photometric data with Sersic extrapolation) and RGGR gravity. The 
black solid line in each of the four large plots is the resulting VD for each model, the yellow dashed 
and blue dotted lines are respectively the stellar Newtonian and non- Newtonian contributions to the 
total VD, and the thin dark gray lines in models iii and iv are the effective error bars for T^ (see text 
for details). The first model (i) assumes isotropy (/3 = 0), the second (ii) assumes /3 € [—1, 1], while 
models iii and iv also consider the expected mass-to-light ratio (T*) in the fitting procedure (see eq. 
(6.2) and Table 1 for details). The vertical dashed line signs the radius above which the observational 
data is considered for the fitting procedure (10 arcsec). The smaller plots at the right show the \a 
and 2a confidence levels for each of the four models. 



The procedure above only evaluates the shape of the VD curve, hence it may lead to 
small values of x^ when the inferred mass-to-light ratios are very far from the theoretical 
expectations. On the other hand, to simply fix the value of T* according to a certain estimate 
seems (at least in many cases) unrealistic, since such estimates have significantly large error 
bars. Here, besides using the standard minimization procedure as above, we also use a novel 
procedure in which the expected value of T* is inserted as an additional data to be fitted, as 
follows. 



N 



x^(T„ei,G2) = j; 



g"modcl(T("*, Ql; 02) ^i) — Cq 



+ N 



Ct 



(6.2) 
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Figure 6. NGC 4374 mass models with MOND. See Fig. 5 for details. 



where T,,, is the theoretically expected stellar mass-to-light ratio, Ct is its error or dispersion, 
and the N multiplying the last term imposes the same weigh for both the curve shape and 
the mass-to-light ratio. 

One way to evaluate the goodness of the fit, which is probably the most straightforward, 
is by the comparison of the x^d values. In the tables 2 and 4 we list the values of the minimum 
of x^ ^iid its corresponding Xrcd- '^^^ i^st is given by the minimum value of x^ divided by 
the difference between N and the number of parameters being fitted. When using the x^ 
given by eq. (6.2), the number of observational data, A^, is added by one.^ 

Besides determining the parameters values for the best fit, we also find the confidence 
level curves for each model. The region inside the n cr's confidence level curve is such that 
x'^ipi, ■■■-.Vk) ^ Xmin + ^X^(^) ^)) where {pi} are the model's free parameters and the values 
of /S.x^{n, k) can be found in standard references on statistics. This is important for testing 
the sensibility of each model and for disclosing correlations between its parameters. 

As a final remark, likewise was done in Refs. [48, 49], the galaxy kinematics inside its 
10" radius is not considered at the fitting procedure. This is done since it is unlikely that 
the Jeans model employed is a reasonable approximation for the galaxy dynamics at a range 



^Another reasonable option would be to add by A*', but either choice would not change our conclusions. 
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Figure 7. NGC 4494 mass models with Newtonian gravity (NG). See Fig. 4 for details. 



SO close to the galaxy center. In all of the VD plots the observational data is present from 
5" , and the fitted curves are accordingly extended. 

6.3 Specific comments on the NGC 4374 and NGC 4494 fits 

NGC 4374 

For Newtonian gravity and isotropic VD (/? = 0), a mass model composed only of stars cannot 
fit the observational data of the galaxy NGC 4374, as it it can be seen in Fig. 4 (see also Ref. 
[48]). Besides being a poor fit considering the shape of the VD curve, it was achieved by using 
a stellar mass-to-light ratio T=k which is considerably above the theoretical expectations. If 
the isotropy condition is dropped, the fit for the shape of the curve can be slightly improved, 
but the anisotropy parameter goes toward the disfavored negative values /? < — 1 and T,,, 
increases even more. This feature is common to other giant elliptical galaxies. 

From Fig. 6 and Table 2, it is clear that MOND fits better the NGC 4374 observational 
data than Newtonian gravity without dark matter. However, it is still a poor fit, in particular 
since: i) There is a significant tendency towards a lower VD curve at large radii (the region 
of PNe data), tendency which is strongly enhanced once the fits considers the expected T,,, 
(the models iii and iv in Fig. 6); ii) if the expected T,,, is not used, the best fit is achieved for 
tangential anisotropy with /3 < — 1. Also, besides these two points, MOND is incompatible 
with the Kroupa IMF expectations for this galaxy, since even for the model iv the derived 
T^, is (far) outside the 2a confidence level (Fig. 6). Other issues of MOND with the giant 
ellipticals can be found for instance in Ref. [53, 65]. 
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Figure 8. NGC 4494 mass models with RGGR. See Fig. 5 for details. 



From Fig. 5 and Table 2, it can be seen that RGGR fit to the data is a satisfactory one 
and outperforms MOND in all the points above. It seems that the single issue that MOND 
does better than RGGR is on it's VD curve continuation towards the galactic center (i.e., the 
extension from 10" (0.83 kpc) to 5" (0.41 kpc)). Since the application of this Jeans modeling 
to a region so close to the galactic center is probably meaningless [48] , that would constitute 
no true advantage to MOND. 

NGC 4494 

There is a wiggle at large radius in our analysis (Fig. 7) that is not presented in Ref. [49]. 
This wiggle is innocuous to the fits, and it appears due to a different convention on the Sersic 
extension starting radius. 

Similarly to other ordinary elliptical galaxies, mass models derived from Newtonian 
gravity and a stellar component provide reasonable fits to the observational VD data [49, 66] , 
se Table 4 and Fig. 7. Actually, the fits with Newtonian gravity without dark matter are 
so good that they have lead to a conflict with CDM expectations [66]. Such lack of dark 
matter- like effects was explained in the context of MOND shortly after [9] , since the typical 
internal accelerations of such galaxies are higher than MOND's ao acceleration. In the case 
of RGGR, considering that NGC 4374 is a much brighter elliptical, the expectation would be 
that NGC 4494's corresponding value of a would be significantly lower than that of the giant 
elliptical, while probably (considering the disk galaxy fits) above O.l/i^. Sharper expectations 
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on the a of NGC 4494 can be drawn after similar galaxies have been analyzed in the RGGR 
framework. The resulting a value in Table 4 lies well inside that region. 

Likewise in the NGC 4374 case, the RGGR derived T* in models i and ii (Fig. 8) is 
naturally compatible with the Kroupa IMF expectations (Table 1). Once such value enters 
as part of the data to be fitted, as in eq. (6.2) and models iii and iv, the parameters 
degeneracies are significantly lowered, and sharper predictions are done. 

7 Conclusions 

In summary, in this paper the RGGR [1, 32] effective mass for stationary spherical systems 
was deduced, we presented general considerations on elliptical galaxies within RGGR, and 
evaluated two specific ellipticals (an ordinary and a giant one) , both with extended observed 
VD by the use of recent planetary nebulae (PNe) data [48, 49]. Good agreement between the 
RGGR VD curve and the observations was found. Considering the galaxies here analyzed, 
no strong tendency towards tangential anisotropy^ was found for RGGR, while this behavior 
appeared in both Newtonian gravity without dark matter and MOND. The stellar mass-to- 
light ratios T* found within RGGR are compatible with the expectations of the Kroupa IMF 
(or other similar IMF that lead to lower T,,, than the Salpeter one), also in accordance with 
the findings of Ref. [1] . 

In the case of NGC 4374, the RGGR fit was clearly better than the MOND's one due 
to a number of features. The discrepancies between MOND and the observational data are 
significantly enhanced when a T^, compatible with the Kroupa IMF is assumed. 

In additional to the standard procedure of using T^, as a free parameter, and only 
evaluating a posteriori whether the fitting procedure yielded reasonable results, additional 
fits were generated in which the expected value of T,,, was part of the data do be fitted. This 
procedure, in comparison with simply fixing the value of T* a priori, is useful to avoid being 
unphysically precise on the value of T* . 

The best way to compare different models is to use the same assumptions and procedures 
whenever possible. This methodology was used here for Newtonian gravity, MOND and 
RGGR. In Refs. [48, 49] the NFW dark matter halo was analyzed for the same galaxies, but 
the higher order Jeans equations together with the kurtosis data were used as part of the data 
to be fitted. In particular, the observational kurtosis can put constraints on the anisotropy 
parameter /?, once the distribution function is assumed to satisfy f{E,L) = fQ{E)L~'^" for 
constant /3 (for further details, see the Appendix B of Ref. [49] and references therein). In 
those references, the NFW halo fits were found to be compatible with the N-body simulations 
expectations and the Kroupa IMF, considering the effects of adiabatic contraction together 
with a crescent anisotropy along the radius (albeit with a low concentration value for NGC 
4494). Comparing the shape of the VD curves, the NFW halo could achieve about the same 
or a better concordance with the observational data than RGGR (remember that NFW 
halos use one more free parameter than RGGR) . It is not easy to draw a straight comparison 
between the anisotropy parameters, but RGGR has shown a tendency towards isotropy while 
being only mildly dependent on the precise value of /3 (since the 1 a uncertainties for RGGR 
are about or larger than 0.5). Had the kurtosis been evaluated as part of the data to be fitted, 
that would probably insert a tendency towards higher radial anisotropy. Finally, considering 
the stellar mass-to-light ratios (T=k), albeit both NFW and RGGR are compatible with the 
Kroupa IMF expectations, RGGR shows a tendency in these galaxies towards lower T,,, 

® Which is usually disfavored, see also Ref. [52]. 
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than NFW. And that may eventually constitute a physical test for RGGR, once the stellar 
populations models become more precise. 

Another well known gravitational theory, f{R) gravity, had its consequences for elliptical 
galaxies recently evaluated in Ref. [67]. The sample of data that we use is similar to 
theirs. The VD curve that they find is similar to the one derived for RGGR, while two clear 
differences are the T,,, of NGC 4374 (RGGR has a lower value, more in the middle of the 
Kroupa IMF expectation) and the /5 of NGC 4494 (since f{R) shows a significative tendency 
towards radial anisotropy). It would be interesting to compare these models using a larger 
sample. 

There is a number of natural developments for RGGR, both from the theoretical side 
and the phenomenological one. In particular, within the context of galaxy kinematics, it is 
desirable to better constrain the a variation. This can be achieved from the analysis of a 
larger sample of galaxies. The a variation with the galaxy parameters is not currently well 
known, but it should be stressed that RGGR without dark matter (or with a relatively small 
amount of it inside galaxies) imply the existence of correlations between dark matter-like 
effects and baryonic matter. The existence of some correlations are well known for a long 
time, in particular here we briefly explored the correlations implied by the fundamental plane. 
It might result that no sufficiently strong correlation compatible with an a as a reasonable 
function of M^:,Re,n to be found, hence galaxy kinematics alone (even if all the individual 
fits of galaxies are satisfactory) may refute, or significantly corroborate, this approach of 
RGGR without dark matter. Other developments at cluster and cosmological scales and on 
lensing effects are being done. 
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A The A variation inside galaxies 

The main purpose of this appendix is to show that, albeit A is not constant inside galaxies 
in the RGGR context, its running is not sufficiently high to lead to any significative changes 
to the galaxy dynamics. The differential equation that governs the A running, eq. (A. 6), is 
below derived from the RGGR action and it is numerically evaluated for the galaxies NGC 
4374 and NEC 4494, Fig. 9. 

With the cosmological 'constant' A, the action (2.3) reads, 

Snoonig] = ^j ^^V^gd'x. (A.l) 

The above, together with some energy- momentum tensor T^^, leads to the following field 
equations (see, e.g. [68]) 

G^, + Ag^, + G DG-i^M. " GV^V.Q-' = ^ T^, . (A.2) 
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From the energy- momentum conservation {T^u = 0), the Bianchi identities {G^u = 0) 
and eq. (A. 2), the divergent of the the terms with derivatives on G in eq. (A. 2) must vanish, 
hence 

V.(^)=^i2V,G-\ (A.3) 

since {OV^ - V^n) G~^ = R^^V^G^^. The above result is exact. 

Equation (A.3) also appeared in Refs. [39, 69]. In [39] the authors use this equation to 
motivate the RG energy scale /i as a function of the Ricci scalar R. The approach we follow 
here is different. Here ^ was chosen from physical considerations in the weak field limit (see 
the comments before and after eq. (2.9)), while the G variation as a function of ^ is given 
by eq. (2.2). 

In order to disclose the variation of A as a function of G and a matter content given by 
a dust component (i.e., (T^) = diag(— pc^ 0)), from the trace of eq.(A.2) we write, 

i?= ^p + 4A + 3GnG-^ (A.4) 

Replacing the above in eq. (A.3), one finds 

V^A= r^^+A + ^GDG-i'jGV^G^^ (A.5) 

The above equation is here interpreted as the differential equation that defines the A variation. 
We do not have an analytical solution for A, but from a given p and from the expression 
of G as a function of the Newtonian potential (i..e, G = Gq/{1 + 2va\n.^N /^q), A can be 
numerically determined. Namely, 



A' = 2.a^ 



— ^ h A + 3ua -— —g- 

c2 V ^N ^% 



(A.6) 



Contrary to G, the derivatives of A do not appear in the field equations (A. 2), hence 
A can be dynamically relevant only if its absolute value is sufficiently high. For comparison 
purposes, the A value as estimated from ACDM cosmology is Aacdm = 1.3 x 10^^^/kpc , 
while AttGop^/c^ ~ (10~^ — 10~^^)/kpc^ inside elliptical galaxies (up to about 5 effective 
radius, with p* the stellar mass, no dark matter). 

To evaluate A inside galaxies in the RGGR context the starting point is to fix the value 
of A "outside" the galaxy, say at 25 effective radius (Re)- To generate the numerical solutions 
presented in Fig. 9, we use that A(25i?e) = Aacdm- One can change the value of A(25i?e) 
by more than an order of magnitude without significative consequences to the numerical 
evaluations here presented. 

In conclusion, within the RGGR framework it is indeed safe to neglect A inside the 
galaxies NGC 4494 and NGC 4374 (and probably the same holds for any other galaxy). 
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